Charge-Carrier Screening in Single-Layer Graphene 
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The effect of cliarge carrier screening on the transport properties of a neutral graphene sheet is 
studied by directly probing its electronic structure. We find that the Fermi velocity, Dirac point 
velocity and overall distortion of the Dirac cone are renormalized due to the screening of electron- 
electron interaction in an unusual way We also observe an increase of the electron mean free path 
due to the screening of charged impurities. These observations help us to understand the basis 
for the transport properties of graphene, as well as the fundamental physics of these interesting 
electron-electron interactions at the Dirac point crossing. 



Unlike normal metals where charge carriers and im- 
purities are highly screened by the Fermi sea, the intro- 
duction of charges to a neutral graphene sheet has sev- 
eral competing effects on its transport properties, due to 
the screening of the electron-electron interaction and the 
screening of long-range impurities (such as charges or va- 
cancies) . While the former is expected to renormalize the 
Fermi velocity and the Dirac point velocity [1| , the latter 
could lead to a decrease of the quasiparticle scattering 
rate or an increase in the electron mean free path. [2, |3| 
These effects are important for those applications, such as 
spintronics, where the electron mean free path is a more 
relevant quantity than the conductivity. [4| Although the 
effect of electronic screening has been extensively studied 
in the past [l|, y, I5l4l5l| , revealing unusual behavior upon 
the introduction of charge carriers on the graphene sheet 
|16| , most of these works have focused on the renormal- 
ization of the Fermi level IT], |l8| , and cannot address the 



important question of how the various electronic effects 
can renormalize the Dirac cone or otherwise contribute 
to conduction. 

Angle-resolved photoemission spectroscopy (ARPES) 
is an ideal tool to probe the electronic properties of 
graphene 19l425[. In many of these works however, the 
starting graphene layer is highly doped, and as a con- 
sequence the effect of charge-carriers on the electronic 
screening is more difficult to discern. Similarly, in the 
case of undoped graphene, the focus has been on how 
the dielectric screening of the substrate has an impact on 
the electronic dispersions 23|, |24[. Therefore, the ques- 
tion of how charge carrier screening affects the Dirac cone 
dispersion and how it differs from dielectric screening re- 
mains an open questions. 

Here we demonstrate the effects of charge-carrier 
screening on a graphene sheet: with the progressive depo- 
sition of small quantities of potassium, we observe a sin- 
gularity in the Fermi velocity and Dirac point velocity; an 
overall renormalization of the valence band; a decrease in 
the quasiparticle scattering rate; and qualitatively differ- 
ent behaviors from the case of dielectric screening. These 
results demonstrate the many ways in which charged im- 
purities can have an impact on the transport properties 



of graphene. 

The h-BN samples were prepared by CVD growth on 
a Cu film, followed by transfer to the h-BN substrate 
and hydrogen annealling. The graphene was placed on 
mechanically exfoliated flakes of h-BN, many layers high 
(opaque to visible light), which were in turn supported 
on a doped Si wafer with native oxide. The sample 
preparation was nearly identical to that described in 
past references (the samples in our experiment were not 
patterned). [H, 0, 29] Our ARPES investigation was 
performed at beamline 12.0.1.1 at the Advanced Light 
Source, at a pressure better than 3xl0~^^ torr, with a 
sample temperature of 15 K, and photon energy 50eV. 
The sample was annealed to 700° C in UHV prior to 
measurement. The sample was electron-doped in situ by 
potassium deposition with an SAES Getters alkali metal 
dispenser at 15K, under which conditions the potassium 
atoms sit above the graphene surface in a disordered ar- 
rangement. [2lJ, |42| The in situ deposition allowed us to 



study the same position on the sample as potassium was 
progressively added. 

Much attention has recently been focused on the 
properties of hexagonal boron nitrid e ( h-BN) as a sub- 
strate for graphene electronics |26l429| . Graphene/h- 
BN has significantly improved transport properties 
and fewer charged impurities than previously stud- 
ied graphene/Si02 systems [30|. ARPES spectra of 
graphene/h-BN are shown in Figure 1. Following the 
maximum intensity, one can clearly observe nearly linear 
energy spectra, characteristic of Dirac electrons |32|. As 
potassium is added to the sample, the Dirac point ap- 
pears and moves to higher binding energy (indicated by 
the black arrows in panels d-h), and the charge density 
of the sample increases. 

As potassium is added the spectral widths do not in- 
crease significantly, which is surprising since impurities 
often broaden ARPES spectra. In Figure 2 we examine 
this effect, showing the widths of momentum distribu- 
tion curves (MDCs, intensity profiles as a function of 
momentum) for different dopings. Panel (a) shows that 
the spectral widths vary almost linearly as a function 
of binding energy, but have an overall offset: the sam- 
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FIG. 1: ARPES spectra of graphene on h-BN. (a) ARPES 
spectrum of single-layer graphene along the F-K direction. 
This is the direction along which all of our data is analyzed 
in subsequent figures, (b) ARPES spectrum along the K-K' 
direction (perpendicular to F-K) shows that the Dirac point 
is at the Fermi level. The spectra in panels (a)-(b) have been 
normalized by the area under the MDCs. (c) The pointlike 
Fermi surface of graphene. (d-h) Doping dependence along 
the F-K direction, with Fermi k-vectors corresponding to kp 
= -0.0035, 0.0037, 0.0126,0.0206, 0.0282 A"\ respectively. 



pie with higher charge density (larger kp) has smaller 
(sharper) MDC widths. 

To quantify this further, in panel (b) we plot the widths 
of MDCs at the Fermi level, as a function of the Fermi 
momentum fcp- The width of an MDC at the Fermi 
level is proportional to the quasiparticle scattering rate 
or inverse mean free path of the photohole[3|, |33|, and 
fcp is proportional to the square root of the charge den- 
sity. We observe that as the carrier concentration in- 
creases in absolute magnitude, the MDC widths decrease. 
This behavior differs from short-range impurity scatter- 
ing, where scattering rates are expected to increase with 
charge density. [3j| The decrease in the MDC widths also 
cannot be due to other many-body effects, which gener- 
ally have vanishing contributions to the spectral widths 
at the Fermi level. [35| Therefore the screening of long- 
range impurity scattering is the only remaining expla- 
nation for the sharpening of spectral features: when the 
increase in charge density improves the screening of long- 
range impurity potentials, the quasiparticle scattering 
rate is expected to decrease with increasing doping. Q 

Figure 2b can be compared to theoretical calculations 
of the scattering rate from long-range impurities [6|, |8|. 
The dashed line shows a fit to a l/fep behavior, where the 
constant of proportionality gives the impurity density: 



ImS ~ a ni,„pUF7r/(2a)/fcF + Const. 
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FIG. 2: The screening of charged impurities can lower the 
quasiparticle scattering rate, (a) The imaginary self-energy 
(ImS), proportional to the width of ARPES spectra and in- 
versely proportional to the quantum lifetime of the photohole, 
is shown for two different dopings (n is proportional to kp). 
In addition to the roughly linear binding energy dependence 
for each spectrum is an overall offset between them: the spec- 
tra with added potassium have sharper spectral features than 
the as-grown sample, a naively counterintuitive finding, (b) 
A comparison of the MDC widths at the Fermi level shows 
that increasing the charge density sharpens spectral features 
by screening the interaction between quasiparticles and impu- 
rities. The dashed line is a l/kp fit to the data, which allows 



us to extract an impurity density of nin 



1.6 X 10' 



Here ImE is the imaginary self energy, a is the effective 
fine-structure constant of graphene, riimp is the impurity 



The inset shows the raw MDCs at the Fermi level, confirming 
that the blue curve is sharper than the red. 



density, vf is the Fermi velocity, /(2a) is a dimension- 
less constant (we use /(2a) ~ 0.22) [1^, and Const is 
an overall offset. Using a = 0.78 (discussed below) and 
vp — 0.85 X lO^m/s (the bare LDA velocity), we find 
"imp = (1-94 ± 0.37) X 10"cm"2, which is typical for 
graphene/h-BN [30|, smaller than the impurity density 
of graphene on Si02 (typically > lO^^cm"^) [36], and an 
order of magnitude smaller than the potassium density 
of the highest doping in figure 1 (1.3 x 10^^, assuming ^ 
1 electron donated per potassium atom [37|). 

To further investigate the effect of long-range screen- 
ing, in Figure 3 we report the doping dependence of 
the graphene bandstructure. The electronic energy- 
momentum dispersion of graphene can be obtained by 
fitting Lorentzian lineshapes to the MDCs and extract- 
ing peak positions as a function of energy. For as-grown 
graphene/h-BN, the extracted dispersion has been dis- 
played as the solid black line in figure 3a, illustrating the 
logarithmic velocity enhancement in the vicinity of the 
Dirac point. 

As the charge density increases, the Fermi velocity de- 
creases. The band dispersions near the Fermi level are 
plotted for several dopings in panel 3b. The Fermi veloc- 
ity, proportional to the slope of the ARPES dispersions, 
decreases by a factor of 4, with a maximum of 2.0 x 10^ 
m/s, as the Fermi k-vector increases by two orders of 
magnitude. These results have been plotted as the red 
open circles in panel 3d and confirm the results of past 
experiments. [17|, llSj 

However, there are several effects that cooperate to 
reduce the Fermi velocity as a function of doping. For 
instance, the bare band of graphene along the F-K di- 
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FIG. 3: Velocity renormalization by charge-carrier screening, (a) An electronic energy-momentum dispersion extracted from 
the peak positions of the MDCs from undoped graphene/h-BN is shown as the solid black curve. The curvature is due to the 
long-range electron-electron interaction [23l ]. and can be compared to the (arbitrary) straight dashed line in the same figure, (b) 
Extracted dispersions show the Fermi velocity as a function of doping. For ease of viewing, the momenta of these dispersions 
have been aligned so that the Fermi k- vectors coincide (the x-axis is given as k-kp). (c) Extracted dispersions show the Dirac 
point velocity as a function of doping. For ease of viewing, the energies of these dispersions have been aligned so that the 
Dirac points coincide (the y-axis is given as E-Ed). (d) The Fermi velocities and Dirac point velocities are given as a function 
of doping (n is proportional to kp). Both the Dirac point and Fermi level show velocity enhancements as the Dirac point 
approaches the Fermi energy. The dashed line is a logarithmic fit to the Dirac point velocities; see text for further details, (e) 
These two Dirac cones give a schematic of where the data is taken in panels (b) and (c). Panel (b) shows dispersions near the 
Fermi level (marked in red in panel e), which change position with respect to the Dirac point as the Fermi level is adjusted. 
Panel (c) shows dispersions near the Dirac point, which remain at the same place on the Dirac cone even as the Fermi energy 
is altered. 



rection is known to decrease in velocity away from the 
Dirac point [38[, with a Van Hove singularity at higher 
dopings. 2l| Given that the band dispersion in panel 



3a demonstrates a logarithmic divergence near the Dirac 
point, even a rigid shift of the Fermi energy would lead to 
a logarithmic Fermi velocity dependence. [23| Finally, the 
electron-phonon interaction is also known to renormalize 
the Dirac cone as a function of doping. [39] 

Therefore, in order to separate the electron-electron 
renormalization from other many-body and bare-band ef- 
fects, in figure 3c we show the electronic dispersions in 
the vicinity of the Dirac point as a function of doping, 
and the velocities have been plotted in panel 3d. [40| The 
Fermi velocity and Dirac point velocity roughly coincide 
for low dopings. However, as the charge density increases 
the effects of the bare band velocity and electron-phonon 
renormalization become significant, causing the Fermi ve- 
locity to be only half as large as the Dirac point velocity. 
Overall, for the purposes of understanding the electron- 
electron interaction, the Dirac point velocity may be the 
most fundamental quantity. [1[ 

The Dirac point velocity can be fitted with a logarith- 
mic dependence [ll, Hsf : 



Wd = "^D 



4(1 



ln( 



Kr 



+ Const. 



(2) 



where Vg = e^ jve gives the value of the dielectric con- 
stant, and the constant a introduces the same fit param- 
eter as Ref I8|, allowing e to effectively increase with 

charge density. From the fit we obtain a — 720A , and 
WDr°/4 = (0.168 ±0.014) x lO^m/s. Using the bare LDA 
value of vd = 0.85 x lO^m/s, we obtain a = 0.78 or 
e° = 3.3. This value of e compares well with the reported 



logarithmic fit to the binding energy dependence in un- 
doped graphenc/h-BN [24j where e° — 4.22; and with the 
expected dielectric screening, being e^ ~ (eh-BN + l)/2 = 
4.02 [41^. 

In figure 4, we observe the differences between two sep- 
arate types of electronic screening effects: screening by 
the graphene charge carriers, which has been modified 
in this study by changing the number of charge carri- 
ers through potassium deposition; and screening by the 
dielectric environment, which can be modified by chang- 
ing the dielectric substrate [24|]. In both cases the Fermi 
level and Dirac point velocities are modified, varying lin- 
early with e, and logarithmically with fep or doping. On 
the other hand, for charge carrier screening the inverse 
screening length qg varies linearly with kp and with 1/e, 
given by %: 



Qs = 4afcp 



(3) 



One might therefore expect differences between these 
screening effects to be observed at high values of mo- 
mentum, where k^qg. 

So in figure 4a we compare the doping-dependence of 
the graphene valence band dispersions over a larger range 
of energy and momentum than figure 3c. At lower val- 
ues of momentum (near the Dirac point), the band ve- 
locities decrease as a function of doping; but at higher 
momenta this trend begins to reverse, with increasing 
velocity as a function of dop ing near 0.1 A^^. In con- 
trast, figure 4b (and Ref. [2J]) shows that when the band 
dispersions of graphene on different dielectric substrates 
are compared, increasing the dielectric constant leads to 
uniformly smaller band velocities at any given value of 
momentum or energy. 
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FIG. 4: Dirac cone renormalization by charge-carrier screen- 
ing, (a) Electronic dispersions near the Dirac point show 
higher binding energy behavior than figure 3c. The inset 
shows the difference between the experimental bands and an 
arbitrary straight line. Panel (b) illustrates the difference be- 
tween charge carrier screening and dielectric screening. Three 
dispersions are shown: (i) as-grown graphene/h-BN; (ii) as- 
grown graphene/SiC(0001); and (iii) doped graphene/h-BN. 
The doping for (iii) was chosen so that the band velocity near 
the Dirac point would match that of (ii). The cartoons in 
panels (c) and (d) illustrate the renormalization effect on the 
Dirac cones for charge-carrier screening and dielectric screen- 
ing, respectively. In panel (c), the renormalization is primar- 
ily restricted to low momenta. In panel (d), the renormaliza- 
tion extends to all momenta within our range of measurement, 
and becoming larger in magnitude at higher momenta. 



To make this comparison more straightforward, fig- 
ure 4b shows three graphene dispersions: (i) as-grown 
graphene/h-BN; (ii) as-grown graphene/SiC(0001); and 
(iii) doped graphene/h-BN. (i) and (ii) have different sub- 
strates and therefore different dielectric environments; (i) 
and (iii) have the same substrate but different charge car- 
rier concentrations; and (ii) and (iii) have different dielec- 
tric screening and charge carrier concentrations. While 
(ii) and (iii) have similar band velocities in the vicin- 
ity of the Dirac point, overlapping for small values of 
k, these dispersions diverge for larger values of k. For 
k^kp the extent of the renormalization (or magnitude 
of the self-energy) is found to be strongly dependent on 
the dielectric constant e, but weakly dependent on the 
screening by charge carriers or kp. This confirms that 
the electron-electron interaction is indeed a long-range 
interaction, with a variable length scale due to the con- 



centration of free charges in graphene. The ways in which 
charge carrier screening and dielectric screening modiiy 
the Dirac cone are illustrated in figures 4c and 4d, re- 
spectively. 

In conclusion, we have demonstrated some of the de- 
tailed ways in which the addition of charge carriers to a 
graphene sheet can have an effect on transport properties 
and the renormalization of the Dirac cone. The electron- 
electron and electron-impurity interactions are found to 
be long-range interactions, and in both cases the addition 
of charge carriers is shown to decrease the length scale 
and strength of the interaction. The increase in charge 
density is also shown to renormalize the Dirac cone in a 
distinct manner from dielectric screening. These results 
illustrate the differences between charge-carrier screen- 
ing and dielectric screening in graphene, illuminating the 
transport behavior of graphene while demonstrating the 
interesting differences between the electronic interactions 
of graphene and those of ordinary metals. 
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